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Large-scale quantum information processors must be able to transport and maintain quantum 
information, and repeatedly perform logical operations. Here we demonstrate a combination of all 
the fundamental elements required to perform scalable quantum computing using qubits stored in 
the internal states of trapped atomic ions. We quantify the repeatability of a multi-qubit operation, 
observing no loss of performance despite qubit transport over macroscopic distances. Key to these 
results is the use of different pairs of 9 Be + hyperfine states for robust qubit storage, readout and 
gates, and simultaneous trapping of 24 Mg + "re-cooling" ions along with the qubit ions. 

The long term goal for experimental quantum information processingis to realize a device involving 
large numbers of qubits and even larger numbers of logical operations [llQ] • These resource requirements 
are defined both by the algorithms themselves, and the need for quantum error-correction, which makes 
use of many physical systems to store each qubit [j], Q . The required components for building such a 
device are robust qubit storage, single and two-qubit logic gates, state initialization, readout, and the 
ability to transfer quantum information between spatially separated locations in the processor [3, Q, E| ■ 
All of these components must be able to be performed repeatedly in order to realize a large scale device. 

One experimental implementation of quantum information processing uses qubits stored in the internal 
states of trapped atomic ions. A universal set of quantum logic gates has been demonstrated using 
laser addressing 0, 0> Q, leading to a number of small-scale demonstrations of quantum information 
protocols including teleportation, dense-coding, and a single round of quantum error-correction [6|. A 
major challenge for this implementation is now to integrate scalable techniques required for large-scale 
processing. 

A possible architecture for a large-scale trapped-ion device involves moving quantum information 
around the processor by moving the ions themselves, where the transport is controlled by time vary- 
ing potentials applied to electrodes in a multiple-zone trap array 0, 0, fioj ] . The processor would consist 
of a large number of processing regions working in parallel, with other regions dedicated to qubit storage 
(memory) . A general prescription for the required operations in a single processing region is the following 
(illustrated in Fig. 1), which includes all of the elements necessary for universal quantum computation 



(1) Two qubit ions are held in separate zones, allowing individual addressing for single qubit gates, 
state readout, or state initialization. (2) The ions are then combined in a single zone, and a two-qubit 
gate is performed. (3) The ions are separated, and one is moved to another region of the trap array. 
(4) A third ion is brought into this processing region from another part of the device. In this work we 
implement in a repeated fashion all of the steps which must be performed in a single processing region 
in order to realize this architecture. 
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to the complete version of record at http://www.sciencemag.org/ The manuscript may not be reproduced or used in 
any manner that does not fall within the fair use provisions of the Copyright Act without the prior, written permission 
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Some elements of this architecture have been demonstrated in previous experiments [g, Il2| . which 
involved transport of ions in a multi-zone trap. However, these experiments did not involve the use 
of techniques required for building a large scale device, limiting the size of algorithms which could be 
performed. Primary limiting factors for these experiments were the loss of qubit coherence, caused by 
interaction with the fluctuating magnetic field environment, and motional excitation, which degrades 
the fidelity of subsequent two-qubit gates because of the finite wavelength of the gate control fields [13 1. 
Motional excitation occurs as a result of imperfect control during transport and noisy electric fields 
emanating from the electrode surfaces [§[. In this work, we store qubits robustly using a pair of energy 
eigenstates in the 9 Be + 2s 2 Si/ 2 hyperfine manifold (shown in Fig. 2) whose energy separation does 
not depend on the magnetic field to first order. For the 9 Be + "qubit" ions used here, this condition 
is met at a magnetic field of 0.011964 T for the "memory" qubit states |1) = \F = l,Mp = 0) and 
|0) = \F = 2, Mp — 1) (The states are labeled using the total angular momentum quantum numbers F 
and Mp). The insensitivity to magnetic field changes is crucial for preserving coherence in the presence 
of ambient temporal field fluctuations [14j . and also greatly suppresses phase shifts caused by spatial 
variations in the average field experienced by an ion as it is transported throughout the multi-zone trap 
array. We remove motional excitation prior to each two-qubit gate by recooling "refrigerant" 24 Mg + ions 
that are trapped along with the qubit ions. Laser cooling this second species sympathetically cools the 
first through the strong Coulomb interaction between the ions ESSE!. 

A benchmark for scalability in this implementation is the repeated performance of a complete set of one 
and two-qubit logic gates combined with quantum information transport. We demonstrate repeatability 
of a unitary transformation U which involves four single qubit gates, a two-qubit gate, and transport 
over 960 /im (the sequence for U is shown in Fig. 3a). Ideally, U implements the operation 
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in the |11), 1 10) , |01), 1 00) basis. We directly compare experimental implementation of U and U 2 using 
quantum process tomography [l9| . Process tomography requires the process under investigation to be ap- 
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FIG. 1: Schematic of the sequence of operations implemented in a single processing region for building up a 
computation in the architecture of 0, A large-scale device would involve many of these processing regions 
performing operations in parallel, along with additional regions for memory. Generalized operations would use 
this block structure repeatedly, with perhaps some of the steps omitted. 



3 



|F=1,M F =-1> 



|1,0) = |1) 



|1.1>=|1 > 



|2,-2> 



12,-1) 



|2,0) 



|2,1) = |0> 



|2,2)=|0 e ) 



FIG. 2: Hybrid qubit storage in the 9 Be + 2s 2 Si/2 hyperfine levels. The states are labeled using the total angular 
momentum quantum numbers F and Mf. |0) are the qubit states used for single qubit gates and transport, 
and |1g), |0g) are used for two-qubit gates. For detection, the |1, —1) and |2,2) states are used. At the applied 
magnetic field (B ~ 0.011964 T), the frequencies for transitions between pairs of states with the same F are well 
resolved. 



plied to sixteen input states, followed by measurement in nine orthogonal bases 2TJ. The input states are 
prepared using a combination of optical pumping and single-qubit operations, with the latter performed 
on each qubit individually. The analysis also requires individual single-qubit rotations, followed by indi- 
vidual state measurement of the qubits. The experiment therefore realizes all of the basic components 
illustrated in Fig. 1. We directly compare U and U 2 by running the experimental sequence for a given 
input/output combination on U and U 2 sequentially (shown in Fig. 3b), making the comparison of the 
two robust against long term drifts in experimental parameters. For each input/analysis combination, 
we repeat this sequence 350 times. 




FIG. 3: a) Schematic of the qubit ion trajectories (solid red and dotted blue lines) and gate operations used 
to implement U. The single qubit rotations are "7r/2" = 7?(7r/2,0) (eq. [2). The two-qubit gate implements 
G = D[(l,i,i, 1)]. b) Full sequence used to perform process tomography on U and U 2 . This sequence is repeated 
350 times for each setting of preparation/analysis. 



The experiment utilizes two 9 Be + and two 24 Mg + ions, trapped in a six-zone linear Paul trap 12|. Each 
9 Be + ion is used to store one qubit, and is accompanied at all times by a 24 Mg + refrigerant ion, which is 
used for sympathetic cooling. The ion order is initialized to 9 Be + - 24 Mg + - 24 Mg + - 9 Be + at the start of 



the experimental sequence, and remains in this order throughout 2l|. 

Coherent manipulations of the internal and motional states of the ions are performed using laser- 
induced stimulated Raman transitions Q. Single qubit gates are implemented in the basis |1), |0) by 
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resonant Rabi flopping, applying the rotation 



3(0/2) -ie-^sin(6»/2)\ 
-ie 1 ^ sin(0/2) cos(0/2) J 

where 9 is proportional to the Raman pulse duration and cf> is chosen by adjusting the relative phase of 
the Raman light fields at the ion. We individually address the two qubit ions by holding them in two 
trap zones 240 /im apart, and switching the laser beams between zones. 

To implement two-qubit gates, we first combine all of the ions into a single zone. The four-ion linear 
chain exhibits four axial vibrational normal modes caused by the Coulomb coupling between ions (2lj |. 
After recombination, these modes contain significant excess energy, mainly caused by imperfect control 
of the potentials used during separation and recombination. Therefore, prior to each two-qubit gate, we 
cool each mode to near the quantum ground state (On) ~ 0.06) using a combination of Doppler cooling 



22 1 . Importantly, the cooling light only interacts 

3- 



and resolved sideband cooling on the 24 Mg + ions 
with 24 Mg + , leaving the qubits stored in 9 Be + intact 

The composite two-qubit gate makes use of a geometric phase gate [7[ to implement G = D[(l, i, i, 1)], 
where D[v\ is a diagonal matrix with the vector v on the diagonal. The phase acquired by the 1 10} and 
1 01} states is obtained by transient simultaneous excitation of the two highest-frequency normal modes 
by use of a state-dependent optical dipole force 0. The state dependence of this force is derived from 
a differential li ght shift between the two qubit states, which is highly suppressed for field-independent 
transitions We thus use a hybrid scheme for qubit storage, mapping the qubits into a different 

state manifold for the two-qubit gate 2^, 24, 25|. Prior to applying the optical dipole force, we transfer 
each qubit into a pair of states with a sizeable differential light shift - the "gate" manifold |1g) = |1, 1), 
|0g) = |2, 2) (Fig. 2). After applying the state-dependent force, we reverse this transfer and the ions 
are again separated [221 ] . The gate manifold is sensitive to magnetic field fluctuations, which can lead to 
qubit dephasing. We suppress these effects using spin-echo techniques [21 1 . 

We employ quantum process tomography to characterize our implementation of the unitary operation 
U, including any experimental imperfections 1^, 2^]. The evolution of the qubit system (including that 



caused by undesired interactions with the environment) is described by a completely positive linear map 



Pout = £jj(pin) 19] on the input density matrix p m — J2i j c i,j I*) 01 > where the Cjj are complex numbers 



and i,j are labels that each run over the eigenstates |11), |10), |01), |00). Following [26j, we represent 
the map by a 16 x 16 matrix 

i% = Ei*>0'i^(i*>oi) ■ (3) 

In order to extract this process matrix, we experimentally apply the process to 16 input states made up 
of tensor products of the states |1), |0), (|0) — |l))/\/2 and (|0) + i |l))/\/2- For each output density 
matrix, we apply nine sets of rotations, which allow us to measure the expectation values of the operators 
a s <8> <Jt, where the <j s ,t run over the Pauli matrices /, a x , a y , <r z . Our state readout performs a projective 
measurement in the Z basis on each ion independently. We first transfer population from |0) to |2, 2), 
and from |1) to |1,— 1), and subsequently drive the cycling transition 2s 2 Si/ 2 |2,2) ^ 2p 2 P 3 / 2 |3,3) 
for 200 /is, where |2,2) strongly fluoresces and 1 1 , — 1} does not [j^ ]. We collect a small fraction of the 
emitted photons on a photomultiplier tube. We run the sequence shown in Fig. 3b 350 times for each of 
the 16 input states and nine measurement rotations. The process matrix is obtained directly from the 
recorded photon counts and measurement /preparation set ting s using a maximum- likelihood method that 
ensures that the reconstructed process matrix is physical [261 ]. 

Experimentally obtained process matrices for one and two applications of U are shown in Fig. 4. 
From the reconstructions, we can calculate various measures of the fidelity with which the processes 
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were implemented. A direct comparison between experimental results and the ideal case is given by the 
entanglement fidelity F = Tr(E ideal E)/16 0- We find Fjy = 0.922(4) for a single application of U, 
and Fjj 2 — 0.853(5) for two applications (error estimates are the standard error on the mean obtained 
from parametric bootstrap resampling [22I]). As an additional measure of operation fidelity, we take the 



mean / of the fidelity /(pidealj Pe) = \^{y/^Jp^^PB^p^S)\ [28( between the output density matrices 
obtained from the ideal and experimental processes for an unbiased set of 36 input states (formed from 
the eigenstates of a s ®a tl where cr S)t run over a x ,a y ,a z ). We obtain a mean state fidelity of = 0.940(4) 
for Ejy and fjy 2 = 0.890(4) for E^ 2 . We can compare these values to the entanglement fidelities using 
the relation / = (4F + l)/5 27|, and see that they are consistent. 




FIG. 4: Reconstructed process matrix for a) U and b) two repetitions of U . The map £{\i) produces a matrix 
£k,i for each element \i) (J\. Hence elements of the matrix E are labeled bym = 4(i- l) + k,n — 4(j — + where 
the factor 4 results from the size of the two-qubit state space. For example, the |11) (00| (i = 1, j = 4) element 
of an input density matrix is mapped to £(|11) (00|), a 4x4 block of E given by m £ [1, 4] and n £ [13, 16]. The 
position of each peak is in agreement with the theoretical prediction. 

To compare the performance of a second application of U relative to the first, we can compare its exper- 
imental repetition £jj 2 (pm) to a numeric repetition of the experimental map £jy(pi n ), i.e., to £jy^(pi n ) = 
£jj{£jj{p- m )). Evaluating the fidelities for each against the ideal case yields F^ 2 /F^jj = 1.003(13) and 
fjj 2 / ffjfj = 1.004(10), indicating that the operation fidelity is the same for each application of U . We can 
also make a direct comparison between the processes performed by our implementation of U and U 2 by 
taking the mean fidelity between = £jy{£jy{p- m )) and p^ 2 = £jy 2 (pin) for the 36 input states. We find 
fiflfjij, Pjjz) — 0.987(3). Although this number is not unity, as might be expected, the deviation can be 
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ascribed to bias in the maximum- likelihood reconstruction method for finite sample size [22j . Our results 
are thus consistent with the same operation being performed by the experiment for each application of 
U. 

Sources of error in our system arise primarily from spontaneous photon scattering (~1.5% per U) [29] ] 
and intensity fluctuations of the Raman light fields at the percent level. In order to characterize the loss 
of fidelity caused by single-qubit rotations, we apply process tomography to the experimental sequence, 
but without the state-dependent force pulses. In this case the ions are always in a product state and the 
process matrix for each can be obtained independently The resulting process matrices have mean state 
fidelity relative to the ideal case of 0.97 for a single run of the sequence (which uses eight rotations per 
ion including qubit manifold transfer and spin-echo pulses). During the two-qubit gate, the spin states 
are entangled with the motion. From separate measurements of motional coherence, we estimate the 
infidelity from this source to be less than 1 x 1CP 3 . 

Many challenges remain before large-scale ion trap quantum information processing becomes a reality, 
including increasing fidelities to those required for fault-tolerant quantum error correction [l|, and 
meeting the considerable technical challenge of controlling ions in large multi-dimensional trap arrays 
flpj ]. Both of these challenges could potentially contain problems which have not been considered here, 
and which may require combining ou r ap proach with alternative methods, for instance entanglement 
distribution using photonic networks [3fJ. Nevertheless, the combination of techniques demonstrated 
here includes all of the basic building blocks required in this architecture, and opens up new possibilities 
for quantum information processing as well as state and process engineering. 
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Bollinger for technical assistance and Y. Colombe for comments on the manuscript. This paper is a 
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Materials and Methods 

The Supporting Online Materials provide further details about the sympathetic cooling, the use of 
state-dependent forces to implement a geometric phase gate using multiple motional modes, transfer 
between the gate and memory qubit manifolds, state detection and error analysis for quantum process 
tomography. 



Sympathetic cooling. 

After recombining the ions in a single trap zone, each of the axial modes of the four ion 9 Be + - 24 Mg + - 
24 Mg + - 9 Be + chain contains significant excess energy. The precise amount is difficult to characterize, since 
the distribution of population in the Fock state basis is not well known. For the experiments presented 
here, we require three stages of cooling to attain the ground state. We start with Doppler cooling, which 
thermalizes the state with a mean vibrational occupation number n ~ 15. We then use cycles of resolved 
sideband cooling on both the second and first sidebands of each motional mode to prepare the ground 
state with a fidelity of ~ 0.94 for each mode. In total, the cooling takes 5.1 ms per U, and this is the 
limiting factor in the operation time (followed by transport/separation, which takes 3.6 ms per U). The 
motional excitation produced in our experiment is much higher than in similar previous work [121 ] , which 
we attribute primarily to use of voltage supplies with an voltage update rate below the typical trap 
frequencies [10( . In future this time should be able to be substantially reduced by using supplies with an 
increased update rate. 
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Two-qubit logic gate. 

The two-qubit gate follows the basic method used in 7] , with the modification that we simultaneously 
excite two motional modes rather than one. This speeds up the gate for a given laser intensity, and 
thus reduces the error due to spontaneous photon scattering. We perform the gate with the ions in the 
spatial configuration 9 Be + - 24 Mg + - 24 Mg + - 9 Bc + . The two modes we excite are the highest frequency 



axial modes, which are separated in frequency by 2n x 251 kHz [2l|. The difference frequency of the two 
Raman light fields that generate the state-dependent force is u) = uj^ + 5 = uj^ — 26, where 103,1^4 are 
the frequencies of the motional modes, and 5 = 2ir x 83.6 kHz. The state-dependent force is applied for 
a duration to = 2tt/5, ensuring that the spin states are disentangled from both motional modes at the 
end of the pulse 7]. We choose the axial confinement such that the two 9 Be + ions are separated by a 
half-integer number of wavelengths of the optical dipole potential. The oscillation of the two 9 Be + ions 
in the LU4 mode is out of phase, resulting in motional excitation of this mode for |0g1g) and |1g0g)- For 
the L03 mode, motion is excited (and thus phase acquired) for all qubit states. This can lead to unwanted 
a z rotations being implemented along with the phase gate. For this reason, and to increase robustness 
against other sources of error (such as Stark shifts) during the motional excitation, we apply the phase 
gate in two pulses, with a R(n, (j)a) pulse on |0g) «-> |1g) applied between the two. This is followed by a 
second R{9, 4>' G ) pulse (see next section). The ratio of contributions to the gate phase due to each mode 
is $ 3 /$4 = 2.72. 



Transfer between the gate and memory qubit manifolds 



In order to implement the state-dependent optical dipole force, we transfer from the field-independent 
manifold |1) , |0) to the "gate" manifold |1g) , |0g). This is performed by applying R(ir, cf)) to |1) |1g) 
and |0) <-> |0g) sequentially (for these and other transitions described in this section, the rotation basis 
is given by (\h), \l}), where h (I) indicates the higher (lower) energy state). Each of these pulses imprints 
a different phase 4>q on the qubit state. After the state-dependent force pulses have been applied, we 
repeat this transfer in reverse order, taking the qubit back into the |1) , |0) manifold and imprinting phase 
4>'i,4>' - While in the gate manifold, we apply two spin-echo pulses R(n, (fie), R(t, </>g) on |1g?) <-> |0g)- 
As a result the final phase imprinted on the qubit by these rotations depends on the phase differences 
<pi — (f>'i, <po — 4>'q and 4>g — <P'g- ^ eacn transition is driven with independent oscillators these phase 
differences are zero. In our case these are non-zero due to the use of tunable frequency sources for driving 
Raman transitions, and we adjust the phase of all subsequent single-qubit gates to compensate for this. 



Detection. 



The state of each ion is read out individually. Due to the lack of multiple detection apparatus, we 
perform state readout sequentially. First, with the ions in separate zones, we transfer population from 
the qubit memory levels into the qubit readout levels |2, 2) and |1, —1), using a sequence of R(tt, 0) pulses 
on |2, 1) — > |2, 2) and |1, 0) — > |1, —1). We then resonantly drive the cycling transition 2s 2 Si/2 |2, 2) <-> 3p 
2 P 3 / 2 |3, 3) for the ion in trap zone A (ion 1) and detect. For an ion in |2, 2), the distribution of the number 
of photons collected in the 200 fj,s detection period closely approximates a Poisson distribution with mean 
10 counts. For an ion in |1, — 1), the distribution has two components, the dominant contribution coming 
from background scatter of the laser light (which has a Poisson distribution with mean 0.4 counts), 
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and a much smaller exponential contribution due to repumping of the |1, —1) state into |2, 2) over the 
detection period (the repumping probability is rj 1 x 10~ 3 over 200 ps). In order to detect the second 
ion, we recombine the ions in zone A, and again drive the cycling transition. To ensure that ion 1 
contributes negligibly to the second detection, prior to recombination we transfer all population in this ion 
to |1, —1) by optical pumping to |2, 2) followed by R(tt, 0) transfer pulses on |2, 2) — * |1, 1), |1, 1) — > 1 1 , 0), 
1 1,0) — > 1 1 , — 1 ) . If the first of these was imperfect, it might leave population in the fluorescing state, 
hence we perform a transfer pulse from |2, 2) — > |2, 1) to protect against this (|2, 1) scatters relatively few 
photons for our detection parameters). 



Tomography error analysis by resampling. 



Statistical errors for values given in the text are derived from parametric bootstrap resampling 31 1. 
Using the process matrix obtained from the experimental data, we generated 100 sample data sets. This 
was performed using a random number generator combined with the known probability distributions 
for the photon counts, and choosing from a normal distribution of Rabi frequencies consistent with 
observed intensity fluctuations for the rotation angles 9 used for preparation/analysis rotations. For 
each sample data set, we obtained a process matrix, and evaluated the relevant fidelities. The quoted 
statistical error estimates are the standard error on the mean obtained from the distributions of results. 
These distributions are centred about one standard error below the value obtained directly from the 
experimental data. This shift is due to the fact that the maximum-likelihood estimation is only an 
unbiased estimator in the limit of infinite sample size, and that the estimation constrains the solution 
to lie within the set of physically allowed process matrices. This implies that the quoted experimental 
fidelities are underestimates. 

The effect of reconstruction bias is illustrated by the following example. We numerically generated 
100 sample data sets using both the experimentally obtained £jy(p) and the application of this same 
map twice £jy(£^(p)). For each data set, we performed the maximum likelihood analysis, obtaining E'^ 
and Ejj 2 - We compared the performance of these two processes by taking the mean fidelity between 
p'yy = £'^(£'^(p)) and p'jy 2 = £'jy 2 (p) over the set of 36 input states. The average value of this fidelity 
over the 100 sample data sets was 0.989(2), where for unbiased reconstruction this value would be one. 
This means that the experimental value of 0.987(3) is consistent with the gate operation being the same 
for each repetition. 



